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Curvature  and  Tangency  Handles  for 
Control  of  Convex  Cubic  Shapes 

G.  Figueroa,  M.  Paluszny,  and  F.  Tovar 


Abstract.  We  consider  the  problem  of  modelling  a plane  convex  shape 
with  a closed  component  of  an  algebraic  non-singular  cubic.  All  nine 
degrees  of  freedom  are  interpreted  as  visual  handles,  namely:  tangency  to 
three  prescribed  lines  at  three  given  points  and  the  curvatures  at  these 
points. 


§1.  Introduction 

Algebraic  curves,  beyond  conics,  where  introduced  in  CAGD  by  Sederberg 
[6]  in  1984.  Paluszny  and  Patterson  [1-3]  studied  splines  constructed  with 
segments  of  cubic  algebraic  curves,  called  A-splines.  Tovar,  Paluszny  and 
Patterson  [4]  looked  at  A-splines  constructed  with  segments  of  singular  al- 
gebraic cubics,  which  are  just  rational  cubics,  with  new,  geometrically  more 
meaningful,  control  handles  for  their  shape.  Non-singular  cubics  are  classi- 
cally well-known  objects  Salmon  [5].  Projectively  they  could  be  of  two  types: 
two  or  one  circuit  cubics.  A two  circuit  cubic  consists  of  two  pieces,  one  of 
which  can  be  realized  affinely  as  a convex  closed  curve,  called  an  oval. 

The  goal  of  this  paper  is  to  study  the  feasibility  of  using  ovals  to  model 
plane  convex  shapes.  We  remark  that  ovals  are  not  splines,  they  are  C°° 
curves.  Therefore,  the  main  advantage  of  using  ovals  for  modelling  of  convex 
shapes  is  that  no  stitching  of  segments  is  required,  which  would  be  the  case 
for  splines.  In  fact,  a convex  shape  represented  by  a cubic  oval  doesn’t  have 
any  joints  at  all.  It  seems  natural  to  place  the  study  of  cubic  ovals  within  the 
context  of  A-splines  because  they  are  connected  components  of  cubic  algebraic 
curves.  Moreover,  the  techniques  to  study  the  shape  handles  for  their  control 
are  similar  to  those  used  for  A-splines  as  mentioned  above.  In  particular, 
we  want  to  control  the  shape  of  an  oval  through  its  control  triangle,  contact 
interpolation,  and  curvatures  at  three  prescribed  points,  see  Figures  1-4. 
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Fig.  1.  Shape  control  by  moving  vertices  of  the  control  triangle. 


Fig.  2.  Shape  control  by  moving  one  of  the  tangency  points. 


Fig.  3.  Shape  control  by  moving  two  of  the  tangency  points. 


Fig.  4.  Sharpening  the  curvature  at  one  contact  point. 


§2.  Barycentric  Coordinates  and  Curvature  at  the  Endpoints 

The  general  algebraic  cubic  in  cartesian  coordinates  x,  y is  given  by 

F(x,  y)  = a30x3  + a03y3  + a2ix2y  + ai2xy2  + a2 oz2  + a02y2 
+ an  xy  + a\ox  + aoiy  + oqo- 
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Fig.  5.  Cubic  with  two  prescribed  tangencies. 


If  the  cubic  interpolates  two  points  ( xo,yo ) and  (22,  Vi),  and  it  is  tangent 
to  two  lines  joining  these  points  to  a third  point  (xify\)  as  shown  in  Figure 
5,  then  its  expression  in  terms  of  the  barycentric  coordinates  (S,  T,  U)  with 
respect  to  the  triangle  of  vertices  (x0,y0),  (xi,yx)  and  (ar2 , 3/2)  reduces  to 

F(S,  T,  U ) = aS2U  + bSU2  - cST 2 - dT2U  + eSTU  + }T3  (1) 

where  a,  b , c,  d,  e and  / are  arbitrary  real  coefficients.  This  was  observed  by 
Sederberg  [6].  In  Paluszny  and  Patterson  [2]  it  was  shown  that  the  curvatures 
ks  and  ku  of  (1)  at  (xo,yo)  and  (X2,  j/2)  are  given  by 

k - c A 

a (V(xi  - xo)2  + (vi  - v o)2)3  ’ 


E _d A 

b (V(xi  - x2)2  + {yi  - i/2)2)3  ’ 


where  A is  the  area  of  the  triangle  with  vertices  (xo,yo),  (%i,  yi)  and  (12,2/2)- 


§3.  Three  Prescribed  Contacts 

We  now  focus  on  the  family  of  cubics  with  three  prescribed  contacts.  Consider 
a triangle  of  vertices  Pq.  Pj  and  P2,  and  let  ( s,t,u ) be  the  corresponding 
barycentric  coordinates.  Choose  one  point  on  each  side,  as  shown  in  Figure  6. 

The  barycentric  coordinates  of  each  Qt  are  Q0  = Qo(0,  t0,  u0),  Q 1 = 
<9i(si,0,«i),  and  Q 2 = <32(52,^2, 0).  The  general  equation  of  a cubic  passing 
through  these  points  and  tangent  (or  singular)  at  them  to  the  sides  of  the 
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§4.  Curvatures  at  the  Contact  Points 

We  need  to  control  the  curvatures  at  the  interpolation  points.  The  first  step 
will  be  to  find  formulas  for  the  curvatures  in  terms  of  the  free  parameters  0.300, 
0030,  «oo3  and  ain  °f  the  cubic,  as  expressed  in  the  barycentric  coordinates 
( s,t,u ).  In  the  next  section  we  will  produce  inversion  formulas  for  the  atjk  in 
terms  of  the  curvatures,  which  will  allow  us  to  find  oval  shapes  with  prescribed 
curvatures  at  the  prescribed  contacts.  It  will  be  convenient  to  express  the 
curvatures  ks  and  ku  at  Q2  and  Qq  respectively,  in  terms  of  the  barycentric 
coordinates  s,t,u  with  respect  to  Po,  Pi  and  P2.  As  remarked  in  the  previous 
section 

1 = C/ 

“ <*/  (P1Q0)3  ’ 

where  A / is  the  area  of  triangle  I , P1Q0  denotes  the  distance  between  Pi  and 
Qo,  and  a/,  bj,  ci  and  dj,  are  the  coefficients  of  the  cubic  in  the  barycentric 
coordinates  with  respect  to  Q 2,  Pi  and  Qo,  compare  with  (1).  And  similarly 

k = dl  A/ 

S bi  (P1Q2)3 

In  fact  if  a,  b,  c,  d,  e,  f are  the  coefficients  of  the  cubic  in  the  barycentric 
coordinates  with  respect  to  the  triangle  P0P1P2  as  expressed  in  (1),  then  the 
coefficients  a/,  6/,  cj  and  d/  of  this  cubic  with  respect  to  the  triangle  I of 
vertices  Q 2,  Pi  and  Qo  are  obtained  by  a linear  change  of  variables: 

_ ito«i(a2  + /32)a 003  — £osi(2a30ot£o,s2si  2toa030^2ui  — 0111*2^2^0^1) 

0/1  *i«l*o 

, _ S2*i(a2  + /32)fl300  — t2Wi(2ao3otosi^2  + 2s2«oaoo3'!ii  — aiiiS2touo*i) 


t\s  lU2 


where  a = 11152^0,  P = wo^2*i,  and 


0030^2  + *2a300 
CI  — J.2  > 


OQ03mQ  + 0030^0 


Using  the  relationship  between  the  areas  of  the  triangles  P0P1P2  and 
Q2P1Q0  as  illustrated  in  Figure  6,  it  follows  that  ku  and  ks  can  be  expressed 
directly  in  terms  of  the  geometry  of  the  triangle  P0P1P2: 

k - 4c/  u°  A (9) 

“ aI  *2  (PqPl)3  ' 


4 dj  S2  A 

V^(ra)3’ 
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Fig.  8.  Curvatures  at  the  contact  points. 

where  A is  the  area  of  the  triangle  PqPiP 2.  To  compute  the  curvature  kt  at 
Qi  with  barycentric  coordinates  (si,0,«i),  we  consider  triangle  II  to  get 

r _ 4d//  t2  A 

It  is  easy  to  show  that  using  triangle  III  instead,  we  obtain  the  same  expression 
for  kt. 

Figure  8 illustrates  the  relationships  between  the  coordinates  of  contact 
points  Qo,  Q 1 and  Q 2 and  curvatures  ks,  kt  and  ku. 

§5.  Inversion  Formulas  and  Shape  Control 

Equations  (2)-(4)  can  be  rewritten  as 

kus\ai  — uqCj  = 0 


kgU^bj  — s$dj  = 0 
ktu\au  — t2cjj  = 0, 

where  the  a,b,c,d  coefficients  are  given  in  terms  of  0300,  ao30>  0003,  and  am 
and  ks,  kt  and  ku  are  proportional  to  the  curvatures  ks , kt  and  ku,  i.e. 


, - l P»P' 

K-u  — Ku  * 


, £ P0P2 

kt  = kt  — T— 


, _ 7 P1P2 

Ke  — Ke  . 


The  formulas  for  the  a; jk  in  terms  of  ks,  kt  and  ku  are 


0300  = u\t\(ktku(a  + /?)3(sis2  - t0u0(a  + /3)ks) 
+ uit2{tlulks  - s%tlku  - slu^kt)), 


«030  = uls\{ksku{a  + (3)3{t0t2  - si«i(a  + (3)kt) 
+ s2u0(ufs3kt  - uftgks  - s3tlku)), 
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aoo3  = tls\{kskt{a  + /?)3(u0mi  - t2s2(a  + j3)ku) 

+ sito(s\t\ku  - tlulks  - slufkt)),  (5) 

am  = (a  + P)3{ksktku(a  + /?)((<*  - (3)2  - 2 a/3) 

4~  ^kgkfUQtQSiU^  4"  2/Cs/q£uoto*2,^2  4“  ^t2S2s-^,u^kfku'j 
- 3(a2  4-  P2)(ht\s\siku  4-  u^t^uiks  + uoS2s\u\kt ) 

+ (a  + f3)((ot  - P)2  + a(3). 

As  established  above,  see  Figure  7,  further  constraints  are  required  for  the 
contacts  to  occur  at  points  of  the  oval,  and  for  the  latter  to  be  contained 
inside  the  triangle  PoPiP2-  Given  the  points  Q0,  Qi  and  Q2,  (2)  determines 
a cubic  which  is  tangent  to  the  sides  of  triangle  PqPiP2  at  these  points.  To 
guarantee  that  the  points  Qo,  Q\  and  Q2  lie  on  the  oval  it  is  sufficient  that  the 
coefficients  <2300,  a.030  and  aoo3  have  the  same  sign.  Indeed,  the  first  condition 
implies  that  the  third  intersection  of  the  cubic  with  each  of  the  lines  s = 0, 
t = 0 and  u = 0 occurs  outside  the  triangle.  Since  each  contact  Q,  accounts  for 
two  intersections,  if  the  cubic  at  every  Qi  bends  inwards  the  triangle  P0P1P2, 
then  the  contacts  actually  occur  at  points  of  the  oval  of  the  cubic  which  then 
has  to  lie  inside  the  triangle.  The  positivity  of  ks,  kt  and  ku  guarantees  that 
the  curve  bends  towards  the  interior  of  the  triangle,  see  [3]. 

Given  the  contacts  at  Qo(0, to,«o)t  <5i(si,0, ui)  and  <32(«2, 0),  see 
Figure  6,  for  0003,  0030  and  aoo3  to  have  the  same  sign,  it  is  enough  to  take 

,0  S2&1 

s t0u0(a  + f3)' 


,0  _ *0*2 

1 uiSi(a f})' 

, 0 _ uQul 

Ku~  s2t2(a+P)' 

Note  that  given  an  oval  with  three  prescribed  curvatures,  when  we  fix 
two  of  them  the  equations  for  the  in  terms  of  the  third  are  linear.  So,  it 
is  easy  to  find  an  interval  range  for  the  modification  of  the  third  curvature, 
while  keeping  the  coefficients  <1300,  0030  and  aoo3  of  the  same  sign. 

§6.  Conclusion 

A convex  shape  can  be  modelled  with  a cubic  oval  controlled  by  a triangle  that 
contains  it.  By  moving  the  vertices,  we  drag  the  oval  along  (see  Figure  1). 
Moreover,  given  any  triangle  and  three  points,  one  on  each  of  its  sides,  we  can 
produce  a convex  shape  that  contacts  the  triangle  at  these  points.  This  convex 
shape  is  our  cubic  oval.  Furthermore,  we  can  prescribe  the  curvatures  at  these 
three  points,  within  precisely  defined  conditions.  If  we  fix  two  curvatures,  the 
third  curvature  can  be  modified  within  an  interval  which  can  be  computed  by 
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solving  a system  of  linear  inequalities.  In  fact  given  the  vertices  of  the  control 
triangle,  the  contact  points  at  its  sides  and  the  three  initial  curvatures  k !?,  kf 
and  k®  as  given  by  (6),  it  is  always  possible  to  modify  any  of  them  keeping 
the  other  two  fixed. 

This  entails  solving  the  linear  inequalities  CI300  > 0,  0030  > 0 and  a 003  > 0, 
or  <3,300  < 0,  <3o30  < 0 and  aoo3  < 0 given  by  the  linear  system  (5). 
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